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Abstract. Using the 15-m Swedish ESO Sub-millimeter Telescope (SEST), CO, HCN, and HCO + observations 
of the galactic star-forming region NGC 6334 FIR II are presented, complemented by [Cl] 3 Pi— 3 Po and 3 P2 — 3 Pi 
data from the Atacama Pathfinder Experiment (APEX 12-m telescope). Embedded in the extended molecular 
cloud and associated with the Hn region NGC6334-D, there is a molecular "void". [Cl] correlates well with 13 CO 
and other molecular lines and shows no rim brightening relative to molecular cloud regions farther off the void. 
While an interpretation in terms of a highly clumped cloud morphology is possible, with photon dominated regions 
(PDRs) reaching deep into the cloud, the data do not provide any direct evidence for a close association of [C i] 
with PDRs. Kinetic temperatures are ~40-50 K in the molecular cloud and <;200 K toward the void. CO and [C i] 
excitation temperatures are similar. A comparison of molecular and atomic fine structure line emission with the 
far infrared and radio continuum as well as the distribution of 2.2 fim H2 emission indicates that the well-evolved 
H 11 region expands into a medium that is homogeneous on pc-scales. If the H2 emission is predominantly shock 
excited, both the expanding ionization front (classified as subsonic, "D-type") and the associated shock front 
farther out (traced by H2) can be identified, observationally confirming for the first time a classical scenario that 
is predicted by evolutionary models of Hn regions. Integrated line intensity ratios of the observed molecules are 
determined, implying a mean C ls O/C lr O abundance ratio of 4.13 ± 0.13 that reflects the ls O/ 17 isotope ratio. 
This ratio is consistent with values determined in nearby clouds. Right at the edge of the void, however, the 
oxygen isotope ratio might be smaller. 

Key words. ISM: abundances — ISM: atoms — ISM: Hn regions — ISM: individual objects: NGC 6334 FIR II - 
ISM: molecules — radio lines: ISM 



1. Introduction 

Young massive stars greatly affect their environment, ei- 
ther impeding or triggering future star formation. Theory 
(e.g., Franco et al. 1990) predicts that soon after the ra- 
diation field of a hot massive star has been turned on, its 

Send offprint requests to: J.S. Zhang, GuangZhou University, 
dennyzhang998@hotmail.com 

* Based on observations with the Atacama Pathfinder 
Experiment telescope (APEX, Chile) and the Swedish-ESO 
Sub-millimeter Telescope (SEST, ESO/La Silla, Chile). 

** Figure IA.1I is only available in electronic form at 
http:/ /www. edpsciences.org 



energetic photons create a supersonic ("R-type") ioniza- 
tion front that moves through the surrounding molecular 
gas, leaving it hot and ionized. The speed of this front 
is continuously reduced by geometric dilution and recom- 
bination. At a time when the speed is not much larger 
than the speed of sound, i.e., when the "initial" Stromgen 
radius is reached, the large pressure gradient across the 
ionization front is supposed to drive a supersonic shock 
wave that sweeps, accelerates, and compresses the sur- 
rounding medium. The ionization front further deceler- 
ates, becomes subsonic ("D-type"), and follows the shock 
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front, while neutral gas accumulates in the interphase be- 
tween the fronts. 

Complications arise when the surrounding gas is 
clumpy. If the ionization or shock front encounters a strong 
negative density gradient, e.g., at the edge of the parent 
cloud, the ionized gas may expand supersonically ("cham- 
pagne phase"). Another complication arises when young 
massive stars form in clusters so that scenarios based on 
a single dominating star may be misleading. Poorly colli- 
mated mainly low velocity stellar winds from young mas- 
sive stars may also contribute significantly to the budget 
of incoming energy (e.g., Tenorio-Tagle 1988; Churchwell 
120021) . 

Since it is known that molecular gas is highly clumped 
(e.g., Wilson & Walmslcy 1989), the less energetic UV ra- 
diation (~6-13.6eV) emitted by the young star(s) may 
leak deeply into the cloud. Such "photon-dominated re- 
gions" (PDRs) may thus make up a large fraction of the 
volume of a star- forming cloud (e.g., Hollcnbach & Ticlcns 
I1997p . The UV radiation from the star as well as associ- 
ated shock fronts heat the dust and the gas. Molecules will 
be photodissociated and surface layers of dust grains will 
evaporate. The particles generated in such ways are ready 
to form new molecular species characteristic of a PDR or 
shocked environment. 

Many star-forming regions have been studied in some 
detail. Usually it is found that the geometry of the scru- 
tinized region is irregular. In the following we present 
a study of an object of particularly simple geometry, 
the periphery of an almost circular region devoid of 
molecular gas. This region forms the eastern part of 
NGC 6334 FIR II, the outstanding source of far infrared 
emission in the prominent southern star forming region 
NGC 6334 (Loughran et al. 1986). Here we present maps of 
the region in a number of molecular transitions (A~3 and 
1.3 mm) that are complemented by profiles from the 492 
and 809 GHz lines of neutral carbon (Ci). Observational 
methods are described in Sect. 2. New data are presented 
in Sect. 3, while Sect. 4 provides physical parameters, in- 
cluding isotope ratios determined from CO. Section 5 sum- 
marizes the main results. 

2. Observations 

2.1. SEST observations 

The J = 1 - lines of 12 CO, 13 CO, C 18 0, C 17 0, HCN, 
and HCO+, as well as the J = 2 - 1 line of 13 CO have 
been observed with the 15-m Swedish-ESO Sub-millimeter 
Telescope (SEST, Booth et al. 1989) in May 1993 (C 18 0) 
and May 1994. For the 3 mm observations, a Schottky 
receiver was used. System temperatures, T sys , were typi- 
cally 480 K on a main-beam brightness temperature (T m b) 
scale. 1.3 mm data were obtained with an SIS receiver with 
T sys ~1600K. 

The chopper wheel method was used for calibration. 
The intensity scale was converted to a T m b scale (see 
Downes 1989; Rohlfs & Wilson 1996) assuming main-beam 



efficiencies of 0.75, 0.70, and 0.60 at v ~ 90 GHz, 110 - 
115 GHz, and 220 GHz, respectively (SEST Handbook, 
1993). The full width to half power (FWHP) beam- width 
was 57" for HCN and HCO+, 45" for the J = 1 - lines 
of CO and its rare isotopomers, and 24" for the J = 2 — 1 
transition of 13 CO. 

The pointing of the telescope was regularly determined 
by mapping the SiO masers toward VX Sgr and AH Sco. 
We estimate the pointing to be correct within 10" (but 
see C 17 in Sect. 4.7), while the calibration uncertainty is 
estimated to be ±15%. As backends, we employed acousto- 
optical spectrometers with 1000 or 1652 contiguous chan- 
nels, a channel spacing of 43 kHz, and a spectral reso- 
lution of 80 kHz (0.27, 0.21, and 0.11 km s" 1 at 90, 112, 
and 220 GHz, respectively). The J = 1 — transitions 
of 12 CO and 13 CO were measured in a dual frequency 
switching mode, with frequency offsets of ±15 MHz. The 
other measurements were carried out in a position switch- 
ing mode, the reference position being at an offset of 
(Aa,A<5) = (—20', 20'). The integration times per point 
were usually 15 to 60 seconds; only for C 17 were in- 
tegration times longer, with typically 5 minutes per po- 
sition. From all the spectra, linear baselines have been 
subtracted. A summary of observational parameters is 
presented in TableQ] Spectra reduced using the GILDAS 
software (e.g., Guilloteau & Lucas 2000) are shown in 
Fig. IA.llF of the appendix that is available electronically. 



2.2. APEX observations 

The 3 Pi— 3 Po and 3 P2— 3 Pi fine-structure lines of neu- 
tral atomic carbon (Ci) were observed in April 2006, us- 
ing FLASH (First Light Apex Submillimeter Heterodyne 
instrument; Heyminck et al. 2006) on the 12-m APEX 
(Atacama Pathfinder Experiment; Gusten et al. I2006P 
telescope. FLASH is a dual-channel heterodyne instru- 
ment operating simultaneously in the 420-500 GHz and 
the 780-880 GHz atmospheric windows. Double side- 
band system temperatures were ~1000-2000K and 3800- 
8000 K on a T m b scale, respectively, mainly depending on 
elevation. The measurements were carried out in a po- 
sition switching mode with the reference position at an 
offset of (Aa, AS) — (20', 0'). On-source integration times 
were 30 to 90 seconds per position; the central position 
was observed repeatedly to check calibration, resulting in 
a total on-source integration time of about 300 seconds. 
Observational parameters are summarized in Table[TJ 

The data were calibrated to a T m b scale adopting main- 
beam efficiencies of 0.60 and 0.43 and forward hemisphere 
efficiencies of 0.95 at v ~ 492 GHz and 809 GHz, respec- 
tively. Full width to half power (FWHP) beam-widths are 
14" and 8" (Gusten et al. 2006). The absolute calibration 



Throughout the paper, B1950.0 coordinates were used to 
facilitate comparisons with previously published results (e.g., 
Brooks & Whiteoak l2"uTjTl Ezoe et al. 120051 and earlier publi- 
cations). 
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Table 1. Observed lines toward NGC6334 



Transition 


Frequency 
(GHz) 


Beam 
(") 


L a) 
(PC) 


b) 


Time c) 
(s) 


SEST 15-m 












HCN(l-O) 


88.632 


57 


0.48 


0.75 


30 


HCO+(1-0) 


89.188 


57 


0.48 


0.75 


30 


12 C 18 O(l-0) 


109.782 


45 


0.38 


0.7 


60 


13 C 16 O(l-0) 


110.201 


45 


0.38 


0.7 


15 


12 C 17 O(l-0) 


112.359 


45 


0.38 


0.7 


300 


12 C 16 O(l-0) 


115.271 


45 


0.38 


0.7 


15 


13 C le O(2-l) 


220.399 


24 


0.20 


0.6 


30 


APEX 12-m 

[Ci] 3 Pi- 3 P 
[C i] 3 P 2 - 3 Pi 


492.160 
809.341 


14 

8 


0.11 
0.07 


0.60 
0.43 


30 
30 



a) To estimate the linear size of the respective beam, 
£>=1.7kpc is adopted (Neckel fWS)) . 

b) Main beam efficiencies; for references see Sect. 2. 

c) On-source integration times: For [Cl], integration times are 
typically 30 seconds per scan, but the central position has 
been observed several times (see Sect. 2.2). 

uncertainty is estimated to be ±25%; the relative calibra- 
tion uncertainty is about ±15%. 

As backend a fast-fourier-transform spectrometer was 
used, containing 16384 channels and covering a bandwidth 
of 1 GHz (Klein et al. 2006). The resulting channel spacing 
is 61kHz (0.04 and 0.02 km s -1 at the lower and higher 
frequencies, respectively). Pointing and focus sources were 
SgrB2(N) and Venus. 

3. Results 

3.1. Molecular line emission 

We mapped the large-scale distribution of the 13 CO J = 

1 — emission over a region of roughly 10'x30' (Fig.QJi) 
covering almost the entire "molecular ridge" of NGC 6334 
that is characterized by a position angle of about 45° (e.g., 
Kraemer & Jackson 1999). At a distance of 1.7(±0.3)kpc 
(Neckel 1978), 1' corresponds to a linear size of 0.5 pc. The 
map spacing is 40" near the ridge and 2' further away from 
the ridge. The northeastern part of the map shows a small 
inconspicious "void" surrounded by extended molecular 
emission. Covering this area that is marked by a square 
in Fig.[lK, maps of 12 CO, 13 CO, C 18 0, C 17 0, HCN, and 
HCO + were taken on a finer grid consisting of 13 x 13 
spectra. These are shown in Figs.QJi-h. The spacing of 
20" implies full sampling for all lines except 13 CO J = 

2 — 1. In these more detailed maps, the void becomes a 
prominent feature common to all molecular lines observed 
by us across the eastern part of the FIR II region. 

Line intensity contrasts are high. Ratios of integrated 
intensity between the position of strongest emission in the 
mapped area (offset relative to the position of minimum 




Fig. 2. Local standard of rest (Vlsr) velocity distribution 
of 13 CO(1-0) across the eastern part of NGC 6334 FIR II 
(B1950.0 coordinates). Contours: -6.5 to -2.0 km s -1 by 
0.5 km s^ 1 . 1' corresponds to 0.5 pc. 



emission: (Acv, AS) — (60", 80")) and the minimum itself 

at a 1950 - 17 h 17 m 26?7, S 1950 = -35° 46' 10" (a 200 o = 
17 h 20 m 47!9; £ 2ooo = _ 35 o 4Q / 08 //) &re 4^ 10; g; 4; 6 ^ and 

6 for CO, 13 CO J=l-0, 13 CO J=2-l, C 18 0, HCN, and 
HCO + , respectively (because of its weakness, C 17 is not 
included here; the error in the ratios can be deduced from 
the 15% calibration uncertainties given in Sect. 2.1 and is 
therefore of order ±20%). The ratios are lower limits be- 
cause not only the main beam but also the sidclobes may 
contribute to the signal seen at the center of the molecular 
void. The size of the region is 1 '.5±0'.2, if measured where 
integrated intensities reach twice the minimum value. This 
holds for CO, for the rare CO isotopomers, as well as 
for the high density tracers HCN and HCO + . The size 
is significantly larger than the beam (Sect. 2.1), so that 
estimates based on 13 CO J=l-0 and 2-1 data (FWHP 
beamwidths: 45" and 24") do agree even without beam 
deconvolution. While in CO the center of the molecular 
void is surrounded by a closed ring of strong CO emis- 
sion, the rare CO isotopomers show an opening of the 
ring toward the south. The high density tracers HCN and 
HCO + also show such an opening toward the south-east 
and east. Here only half of the ring is prominent, even 
though the minimum is clearly identified. From this we 
conclude that molecular column densities are particularly 
high in the east, north, and west, while densities reach the 
highest values near the northern and western edges of the 
molecular void. 

Calculating the mean velocity (first moment) of the 
observed emission, we show the velocity distribution of 
the 13 CO J=l-0 line in Fig.d The other lines show 
a similar distribution. Molecular gas with "low" veloci- 
ties dominates the northeast (down to -6.5kms _1 ), while 
gas with "high" velocities is observed in the southwest 
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Fig. 1. Contour maps of integrated intensity, J Tmb<1v. V corresponds to 0.5 pc. The measured positions are in- 
dicated by markers, a) Large-scale 13 CO J=l— map. The lowest contour is 8Kkms _1 , the contour increment is 
18Kkms _1 , and the rms noise is 1.11 K km s -1 . The arrow points toward the area shown in the other maps of this 
figure, including the emission of seven molecular lines centered at 011950 = 17 h 17 m 26!7, <5i95o = — 35°46'50" (02000 = 
17 h 20 m 47!9, <5i95o = -35°49'48"). b) 12 CO J = 1 - (lowest contour: gOKkms" 1 , increment: 33Kkms-\ rms 
noise: 5.85Kkms _1 ). c) 13 CO J =1-0 (lowest contour: 8Kkms _1 , increment: 18Kkms _1 , rms noise: 1.82Kkms _1 ). 
d) 13 CO J=2-l (lowest contour: SSKkms" 1 , increment: 21Kkms" 1 , rms noise: 0.91 Kkms" 1 ). e) C 18 J=l-0 
(lowest contour: 0.5Kkms _1 , increment: 2Kkms _1 , rms noise: 0.52 K kms _1 V f) C 17 J=l (lowest contour: 
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Fig. 3. Spectra of the [C i] 3 Pi- 3 P 492 GHz line (chan- 
nel spacings 0.15kms _1 (left panel) and 0.30kms _1 (right 
panel)) and the 3 P 2 — 3 Pi 809 GHz line (channel spacings 
0.68kms _1 (left panel) and 0.90kms _1 (right panel)) to- 
ward the molecular void region. In the left panel, west 
is up and east is down; in the right one north is up 
and south is down. The central position is at CH1950 = 
17 h 17 m 26?7, £1950 = -35°46'10" (a 2000 = 17 h 20 m 47!9, 
<5 2000 = -35°49'08"), the spacing is 10". The ordinate is in 
units of main-beam brightness temperature (K). To facil- 
itate comparisons, Vlsr = -5kms _1 is marked by dashed 
vertical lines. Adopted rest frequencies are 492.160651 and 
809.341970 GHz. 



(up to -2 km s 1 ) . The velocity gradient is approximately 
5kms _1 pc _1 ; linewidths are ^5-6 kms . 



3.2. Atomic carbon emission 

To study the neutral atomic component of the molecular 
void, we observed the two [C 1] lines of the 3 P fine structure 
system at 492 ( 3 Pi- 3 P ) and 809 GHz ( 3 P 2 - 3 Pi). [Ci] is 
considered to be a classical PDR tracer (e.g., Hollenbach 
& Tielens H997p . The cross consisting of 19 positions along 
the E-W and N-S directions is shown in Fig.[3J The spacing 
is 10". While one might have expected particularly strong 
emission at the inner portion of the molecular shell, with 
weaker emission deeper inside the molecular cloud, this 
is not seen. Instead, the two atomic carbon lines show 
a trend that is similar to that of the molecular emission 
lines. The line intensity increases when leaving the inner 
region of the void, reaching a high plateau where the sur- 
rounding molecular cloud is located. 

As in the case of the rare CO isotopomers, the molec- 
ular edge toward the south exhibits weaker emission than 
the positions at the other ends of the cross. Because of the 
much smaller beamsize of the [C 1] measurements when be- 
ing compared with the molecular data (Sect. 3.1), the size 
of the region devoid of strong emission as well as spa- 
tial fine structure can be evaluated with higher accuracy. 
Along the E-W axis, we find an extent of 80" for the molec- 
ular void; along the N-S axis the corresponding value is 
100", defined in the same way as for the molecular emis- 
sion in Sect. 3.1. Toward the west, the edge of the void is 
particularly pronounced. Here the [C 1] emission rises by a 
factor of > 2 within 10" (0.1 pc). This is the region that 
appears to be characterized by high column densities and 
high particle densities (Sect. 3.1). A significant difference 
between the distributions of the 492 and 809 GHz lines of 
[Ci] is not apparent. Within the limits of spatial resolu- 
tion, the morphology of [C 1] is compatible with that of 
the molecular tracers. 

For the central position, the main beam brightness 
temperatures are ~4 and ~6K for [Ci] 3 Pi— 3 Po and 
3 P 2 — 3 Pi, respectively, while temperatures of order 20- 
30 K are reached in the outermost eastern and western 
positions. The contrast between the molecular edge and 
the central emission features is ~5±1. This is consistent 
with the corresponding molecular ratios choosing the same 
reference positions. While the mean velocities agree with 
those of the molecular lines (see Fig. [2]), the [Ci] lines 
are slightly narrower, typically ~4kms _1 . This is likely 
caused by the different beamwidths (Table [T]). 



4. Analysis and discussion 

In the following we give a brief overview on NGC 6334 with 
an emphasis on FIR II. Then physical parameters of the 
gas near the void are evaluated, starting with an analysis 
of [C 1] that is followed by CO and the high density tracers 
HCN and HCO + . After searching for traces of molecular 
outflows seen in CO, 18 0/ 17 ratios and the evolutionary 
stage of the entire region are also discussed. 
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4.1. NGC6334: A brief overview 

NGC 63341 hosts some of the most active sites of star for- 
mation in the Galaxy. Embedded in an extended region 
of Ha emission that encompasses almost 1° (Meaburn & 
White 1982; Burton et al. I2000[) lies a giant molecular 
cloud (GMC). Its ridge of length ~30' at P.A. - 45° (see 
Figs.QJi and[4ji for 13 CO J— 1-0 emission) shows spectac- 
ular, complex structure including dense molecular clumps, 
maser sources, Hn regions, and luminous sources of dust 
emission (Fig.|4|). The ridge has been detected in many en- 
ergy bands ranging from low frequency radio waves up to 
80 keV (Burton et al. l20001 Caproni et al. [20001 Kraemer et 
al. [20001 McCutcheon et al. [20001 Persi et al. [20001 Sandell 
[20001 Sarma et al. [20001 Brooks & Whiteoak[20QU Carral 
et al. [20021 Bykov et al. [20051 Ezoe et al. [20061 Hunter 
et al. IMJSl Leurini et al. [217051 Schilke et al. [217051 and 
references therein). 

Along the active ridge, NGC 6334 FIR II is the 
outstanding source of far infrared emission (Lfir ~ 
3xl0 5 L Q ; see Fig.0t). Its peak shifts with increasing 
wavelength to the west, with the 21/im and 134/im peaks 
being displaced by about 2' (Loughran et al. [1986). 
Slightly southeast of the infrared peak lies an H II region 
(NGC 6334-D) of diameter 2' that is known to show an al- 
most circular, featureless thermal radio continuum distri- 
bution (Brooks & Whiteoak 2001 and references therein) 
and that is characterized by an exceptionally low extinc- 
tion (the "molecular void"; see Straw & H viand fl989al for 
an extinction map). This Hn region appears to be pow- 
ered by a highly obscured 06.5 star, possibly the most 
massive star of the entire GMC fNeckel [19781 Rodriguez 
et al. [T982l Loughran et al. [T986l Straw et al. H989cj) . 

Straw & Hyland (|1989b|) imaged the H 2 2.12/zm v=l-0 
S(l) line that shows an intense arc at the southern, east- 
ern, and northern edges of the radio continuum source 
(see also Fig. [4b ) being part of a ring surrounding the Hn 
region. In addition, they took near infrared broad band 
spectra toward selected positions. Analyzing line ratios 
they conclude that most of the H 2 emission is caused by 
shocks rather than by fluorescence triggered by UV radi- 
ation. However, extinction is severely affecting measured 
line ratios, so that the case is not entirely clear (see also 
Burton et al. 2000). The good correlation between this H 2 
ring and the dust extinction led Straw & Hyland (1989b) 
to propose that young stars have swept the inner part of 
the ring (the molecular void) free of gas and dust and that 
the H 2 emission is powered by the interaction between the 
stellar wind and the remaining molecular gas. The shell is 
also seen in the 3.3/im emission of polycyclic aromatic 
hydrocarbon (PAH) (Burton et al. l2"000j) . Overall, FIR II 
appears to belong to the more evolved centers of activity 
in NGC 6334 because the far infrared source and the as- 

2 The nomenclature of molecular, maser, and infrared fea- 
tures is quite complex (see the appendix in Kuiper et al. 1995) . 
We follow recent publications in using the nomenclature given 
by McBreen et al. (119790 and Gezari (119820 for the infrared 
sources and Rodriguez et al. (1982) for the Hn regions. 



sociated H II region are more extended than other similar 
targets in the region (Loughran et al. 119861 Straw et al. 
1989c) and because there is no indication of particularly 
low molecular densities in the surrounding environment 
(e.g., Sect. 3.1). 

4.2. Global distribution of the molecular gas 

The distribution of our large-scale 13 CO J=l-0 map 
(Fig.[T]i) agrees well with the CO J=2-l and CS J=3-2 
maps presented by Kraemer & Jackson (j!999|) . Our 13 CO 
J=l-0 map shows pronounced peaks toward NGC 6334 
FIR I (up to -ISOKkms- 1 ), IV and V, while FIR II and 
III are associated with tongues of ridge emission extending 
southwards from FIR I and northeastwards from FIR IV, 
respectively. The less well sampled southwestern part of 
our 13 CO map, not covered by Kraemer & Jackson (1999), 
shows no significant additional features. For the narrow 
ridge, our map is also consistent with the 13 CO J=2-l 
map of Kraemer & Jackson (1999; their Fig. 3) that cov- 
ers a smaller region. 

4.3. [Cl] excitation, optical depth, column density, and 
correlation with 13 CO 

The C I triplet with lines near 492 and 809 GHz is a useful 
tracer of excitation. Our measured [Cl] 3 P 2 — 3 Pi bright- 
ness temperatures are slightly higher than those of the [C i] 
3 Pi— 3 P line (see below). This is consistent with results 
from other galactic star- forming regions (e.g., Schneider et 
al. 2003). Obtaining Gaussian fits to the [Cl] spectra, we 
calculated integrated line intensities for five key positions, 
the central one and the four outermost positions along the 
four cardinal directions (see Fig. [3]). Assuming local ther- 
modynamic equilibrium (LTE), optically thin line emis- 
sion, and uniform cloud coverage, the integrated intensity 
ratios allow us to determine the excitation temperature 
and optical depths of both transitions as well as the beam 
averaged C I column density (e.g., Stutzki et al. H997j Weifi 
et al. I2003[) . Table [2] displays the resulting parameters and 
the equations used. The ratios of integrated line intensities 
(~1.2-1.5) and excitation temperatures (50-70 K) are sim- 
ilar for the four edge positions. Toward the central position 
of the molecular void, however, the line ratio is almost two 
and the excitation temperature must be <^200 K. This is 
in accordance with the temperature of the photodissoci- 
ated gas assumed by Kraemer et al. (|2000p . 250 K. While 
in view of possible systematic calibration errors the ex- 
act value given in Table [2] may be debatable, we note 
that the integration time toward the central position is 
much higher than toward all other positions, so that a low 
signal-to-noise ratio does not affect the deduced parame- 
ters. Furthermore, all adjacent positions still belonging to 
the inner part of the void confirm, albeit with lower signal- 
to-noise ratios, the high [Cl] 3 P 2 - 1 P / 3 Pi- 3 Po line in- 
tensity ratio. Since the [C i] lines are likely thermalized 
(T ox ~ Tkin), we conclude that the inner part of the stud- 
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Table 2. [Ci] excitation parameters for five key positions (center and outermost positions toward E, S, W, and N; see 
Fig.©* 



Pos. 


f T mb dv 




7? 


7cx 


T 1-0 


T2-1 


Nci 




(1-0) 


(2-1) 


(1-0) 


(2-1) 
















(K) 








(K) 






(10 17 cm- 2 ) 


(0,0) 


14.3±1.7 


28.2±4.8 


4.1±1.0 


6.5±1.7 


1.97±0.30 


>200 


<0.02 


<0.04 


~2 


E 


80.9±3.9 


79.2±13.5 


18.6±4.7 


27.1±6.8 


0.98±0.15 


51±13 


0.6(+0.5,-0.2) 


£0.8 


~15 


W 


70.3±3.7 


72.1±13.6 


18.1±4.5 


26.6±6.7 


1.02±0.15 


53±9 


0.6(+0.2,-0.1) 


1.3(+2.2,-0.4) 


~13 


N 


92.4±4.6 


113.3±13.9 


16.5±4.5 


21.0±5.3 


1.23±0.18 


72±20 


0.3(+0.2,-0.1) 


0.5(+0.4,-0.2) 


~15 


S 


40.6±4.5 


39.2±14.0 


11.0±2.8 


13.7±3.4 


0.96±0.14 


49±10 


0.3(+0.2,-0.1) 


0.6(+0.3,-0.2) 


~6 



* Columns: 

Column 1: Offset position; E, W, N, and S denote the outermost positions shown in Fig.[3j with offsets of 90" along the four 
cardinal directions. 

Columns 2 and 3: Integrated main beam brightness temperatures of the 3 Pi— 3 Po and 3 P2— 3 Pi lines; errors are determined 
from Gaussian fits. 

Columns 4 and 5: Main beam brightness temperatures of the 3 Pi— 3 Po and 3 P2— 3 Pi lines; errors reflect uncertainties in absolute 
calibration (±25%; see Sect. 2) and are not providing information on signal-to-noise ratios. 

Column 6: Integrated intensity ratio ([Ci] 3 P2 — 3 Pi / [Ci] 3 Pi— 3 Po) with 15% uncertainties as estimated in Sect. 2. 
Column 7: Excitation temperature, derived from T cx = 38.8/ln(2.11/i?) (Schneider et al. 2003; for R, see Col. 6). 
Columns 8 and 9: Optical depths, derived from ri-o = -ln(l-T mb (l-0) x (exp(23.62/T cx )-l)/23.62) and r 2 _i = -ln(l-T mb (2-l) 
x (exp(38.8/T cx )-l)/38.8) (Schneider et al. 2003). 

Column 10: [Ci] column density, derived from Nci = 5.94xl0 15 x ri_ /(l-exp(-ri_ )) x [(1+3 exp(-23.6/T cx ) + 5exp(- 
62.4/T cx ))/3exp(-23.6/T cx )] x fT mb (l - 0)dv (see Schneider et al. 2003). 



ied region is dominated by gas with temperatures that are 
much higher than those in the predominantly molecular 
outer regions. While the [C i] optical depth at the center 
of the void is one to two orders of magnitude below those 
at the molecular edges, the total C I column density is as a 
consequence of the high excitation, only a factor of several 
smaller. 



With respect to intensity and morphology, the 492 and 
809 GHz lines of neutral carbon are known to be corre- 
lated with low- J CO emission (e.g., Ikeda et al. I2002p . 
The [Ci] 3 Pi- 3 P and 3 P 2 - 3 Pi fine structure lines have 
upper energy levels of 23.6 K and 62.5 K above the ground 
state, respectively. Critical densities are of order 10 3 and 
10 4 cm~ 3 for collisions with H2. These critical densities 
agree with those of the lowest rotational transitions of 
optically thin CO isotopomers (e.g., Mao et al. [2000). 
Toward the molecular void, main beam brightness temper- 
atures of the [C 1] 3 Pi— 3 Po line are similar to those of the 
13 CO J=l-0 transition and about a factor of ~1.6 lower 
than those of 13 CO J =2-1. Correlations between 13 CO 
and [C 1] emission are shown in Fig.[5l Ignoring differences 
in beam size we compare the integrated intensities for 
the 17 positions with both 13 CO and [Ci] emission mea- 
sured. Results from linear fitting are J(CI( 3 Pi — 3 Po)) = 
(0.89 ± 0.11)/( 13 CO(J = 1-0)) + (-1.00 ± 4.67) and 
7(CI( 3 Pi - 3 P )) = (0.47 ± 0.06)/( 13 CO(J = 2-1)) + 
(—9.12 ± 6.10), with correlation coefficients of r=0.88 and 
r=0.90, respectively. Again, this is consistent with results 
from other molecular clouds (e.g., Plume et al. 119991 Ikeda 
et al. 12002)) . 




(K km/s) 



l( CO(J=2-1)) 



Fig. 5. 

the [Ci] 



A comparison of integrated intensities between 
3 Pi- 3 P and the 13 CO J=l-0 and 2-1 lines. 
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Fig. 6. Integrated 12 CO/ 13 CO J=l-0 line intensity ratios 
near N6334-D. Contours: 4-9, increment: 1. 



-35°45'00" - 



35°46'00" 4 




17 17 30' 



Fig. 7. Integrated 13 CO J=2-l/J=l-0 line intensity ra- 
tios near NGC6334-D. Contours: 1.5-3.0, increment: 0.3. 



4.4. CO 

Not only atomic but also molecular lines can be a measure 
of excitation conditions. Figures [H] and [7] show the inte- 
grated CO/ 13 CO J=l-0 and 13 CO J=2-l/J=l-0 line ra- 
tios for the area covering the molecular void. 12 CO/ 13 CO 
ratios peak, not surprisingly, at the very center of the map 
where the expanding H n region has removed most of the 
molecular gas along the line of sight. Here the line ratios 
almost reach a value of 10, while values of —4 character- 
ize the partial ring of relatively high column density that 
opens up toward the south (see also Sect. 3.1). 

CO must be optically thick because CO/ 13 CO line in- 
tensity ratios are much smaller than the "local" 12 C/ 13 C 
ratio of —70 (e.g., Wilson & Rood 11994]) . From the peak 
CO temperature, assuming uniform coverage of the beam, 



the CO excitation temperature T . 

by 



can then be estimated 



hi/ 



mb.CO 



(e h " 



/kT OJ 



1 



(e h " 



/kT c , 



1 



With T mb - 40-50 K (T cmb =2.73K) in the molecular 
cloud surrounding the void we get T cx values that are of 
the same order as T m b. CO must be thermalized because 
high optical depths reduce the low critical densities even 
below the values mentioned in Sect. 4.3. Being thermalized 
and optically thick, the measured T^co values provide 
an estimate of the kinetic temperature of the molecular 
gas. Tkin ~ 40-50 K. In the case of self-absorption, we 
would have only obtained lower limits to T cx and Tkin- 
The CO profiles (Fig. lA.l]) show, however, no evidence for 
this effect. 

In the Ci analysis presented in Sect. 4.3 (the E, W, N, 
and S positions of Table [2|) we obtained similar excitation 
temperatures as for CO. Does the good agreement be- 
tween C I and CO excitation provide a further argument 
that CO self-absorption can be neglected, that C I may 
also be thermalized, and that the APEX and SEST beams 
are approximately uniformly filled with C I and molecular 
gas, respectively? 

Assuming that (1) the 13 CO lines are optically thin, 
that (2) the cloud covers the beam uniformly, and (3) 
neglecting the microwave background (a good approxi- 
mation for Tex^lOK), a comparison of the 13 CO J=2-l 
and 1-0 lines provides yet another measure of excitation. 
Here again line intensity ratios peak near the center of 
the molecular void, reaching ~4, while the ratio decreases 
to —2 in the vicinity of the partial ring of high column 
density (Fig.[7|). We then obtain with 



T, 



ex,13CO 



10. 



hi 



T, 



nib, 2-1 



4 x T mb i. 



excitation temperatures near infinity at the center of the 
void (within the errors consistent with the high excitation 
obtained from C i; see the (0,0) position in Table [2J and 
values — 15 K near the "molecular ring" . These latter tem- 
peratures disagree with those obtained from CO and Ci. 
Is this due to subthermal excitation of 13 CO, is it caused 
by the fact that Ci traces mostly the UV irradiated sur- 
faces of clouds, or is it an optical depth effect? If 13 CO 
is not optically thin, our derived value for T eXi i3co would 
only be a lower limit. 

Since the critical densities of CO, 13 CO, and Ci are 
similar, different excitation temperatures must either be 
related to different spatial distributions or to optical depth 
effects. In the case of line saturation, photon trapping 
slows radiative rates by — r _1 , leading to a greater dom- 
inance of collisional processes that tend to shift excita- 
tion temperatures upwards, toward Tki n - With a CO/ 13 CO 
abundance ratio of 70 and line intensity ratios of order 4, 
r(CO J=l-0) -20 and r( 13 CO J =1-0) -0.3. Applying an 
LVG code, with H2-CO collision rates from Flower (|200ip 
and an ortho/para H2 abundance ratio of 3 (not a critical 
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parameter), we find a density of n(H 2 ) ~10 3 cm 3 and an 
optical depth of order 1.5-2.0 for the 13 CO J=2-l line. 
Thus the 13 CO J =2-1 line is likely moderately optically 
thick and the equation for T ex> i3co indeed provides too 
low a value. In addition, 13 CO is not fully thermalized 

(T ex ,13CO < 7c Xi i2Co)- 

To summarize, there is evidence that C I and CO show 
excitation temperatures of T ex ~ 40-50 K, which is com- 
parable to the kinetic temperature. iV( 13 CO) ~ 10 17 cm~ 2 
and n(Ha) ~ 10 3 cm~ 3 . These values hold, however, only 
for the parent molecular cloud. Within the molecular 
void, excitation and kinetic temperatures are far higher 
but too uncertain to derive reliable (column) densities. 
According to Straw & Hyland (|1989al their Fig. 9), A v 
<10 m in the void, corresponding to a column density of 
N{R 2 ) < 3xl0 22 cm- 2 (Bohlin et al. 11975)1 or iV( 13 CO) < 
4xl0 16 cm- 2 ([H 2 ]/[CO] = 10 4 , [CO]/[ 13 CO] = 70). IRS- 
24, likely the star that is responsible for the molecular 
void, is reported to show an extinction of <;28 m (Straw 
et al. 11989c")) . corresponding to a 13 CO column density of 
N( 13 CO) ^10 17 cm~ 2 . This is compatible with the char- 
acteristic value already calculated for the molecular cloud. 
Apparently, the star is situated near the western edge of 
the void where column densities are high (IRS-23, a little 
nearer to the void, shows A v ~ 15 m ; Straw et al. 1 1989c]) . In 
view of the high obscuration of IRS-24, the star appears 
to be located behind the bulk of the molecular gas. 

4.5. Other molecules 

The excitation temperatures calculated for CO and its iso- 
topomers cannot be used for HCN and HCO + , which have 
much higher critical densities. An analysis of the HCN 
hyperfine components shows that the HCN J =1-0 transi- 
tion, wherever it is strong enough to show its characteristic 
triple, is optically thin. Adopting T cx = 10 K (a reason- 
able guess implying subthermal excitation), characteris- 
tic column densities are of order JV(HCN) ~ 7xl0 13 cm -2 
and iV(HCO+) ~2xl0 13 c m - 2 for the molecular cloud. To 
reach T m b values of 5 K (or more in a few positions), LVG 
calculations with Tkin = 45 K and collision rates of Schoier 
et al. (|2005|) require n(H 2 ) k, 10 5 cm" 3 . This is a much 
higher density than that derived from CO, indicating a 
high degree of small-scale clumping (see also Sect. 4.8). 
The HCN/HCO+ J=l-0 line intensity ratio is displayed 
in Fig.[8j In the molecular cloud, HCN is stronger, in par- 
ticular in the south-western region that forms a part of 
the molecular ridge. Signal-to-noise ratios are low for both 
lines near the void, so that their ratio shows a large scatter. 
Nevertheless, an inspection of the individual profiles con- 
firms the impression given by Fig.[8j that HCO + becomes 
about as strong as HCN. This may be a consequence of rel- 
atively high kinetic temperatures, strong radiation fields, 
and low densities near the edge of the molecular cloud 
(e.g., Fuente et al. [15551 Chin et al. fTMTl Brouillet et al. 
120051 Christopher et al. 12005)) . 
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Fig. 8. Integrated HCN/HCO+ J=l-0 line intensity ra- 
tios observed with a beam size of 57" toward NGC 6334-D. 
Contours: 0.7-2.5, increment: 0.2. 

Table 3. Average integrated line ratios toward NGC6334- 
D a ) 



Transitions Ratio 

/(CO J=1-0)//( 13 CO J=l-0) 5.3±0.1 

J( 13 CO J=2-l)//( 13 CO J=l-0) 2.2±0.1 

/(CO J=l-0)//(HCN J=l-0) 7.1±0.2 

/(HCN J=1-0)//(HCO + J =1-0) 1.7+0.1 



a) Given errors were obtained from Gaussian fits (standard 
deviation) . 



Averaged integrated intensity ratios are summarized 
in Table [3) The intensity ratio of HCN to CO could be 
a qualitative and even a quantitative measure of pressure 
(Heifer & Blitz ll99"7j) . Our mean values are consistent with 
observational results from other clouds (e.g., Pirogov 1999; 
Heifer & Blitz HMJ) . Unlike the HCN and HCO+ J=l- 
lines, CO J =1-0 is usually opaque tracing lower den- 
sity gas. Hence a stronger correlation is expected between 
HCN and HCO+ than between CO and HCN or HCO+. 
This is confirmed by Fig.[9j 

4.6. Outflows in NGC 6334 FIR II 

Straw et al. (1989b) found an arc of shocked vibra- 
tionally excited H 2 enveloping the H II region NGC 6334-D 
(Sect. 4.1) that coincides with the molecular void. This 
was interpreted in terms of molecular gas swept away 
by an outflow powered by the central star(s). Brooks & 
Whiteoak (|2001[) mention that so far no further obser- 
vations have been made to confirm the presence of this 
outflow. Although we have not tested this idea by directly 
observing shocked H 2 with high velocity resolution, we 
have collected well- resolved spectra from the CO J =1-0 
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NGC6334 CO(l-O) 




10 20 30 40 50 60 70 80 90 

l(HCN) (in K km/s) 




l(HCN) (in K km/s) 



Fig. 9. Comparison of integrated intensites of the J=l- 
transitions of CO, HCO+, and HCN. The lines mark 
least square fits with /(CO J=l-0) = (3.55±0.24) J(HCN) 
+ (123±10) (correlation coefficient: 0.75) and /(HCO+ 
J=l-0) = (0.56±0.02) /(HCN) + (1.47±0.89) (correlation 
coefficient: 0.90). 



transition, which is the outflow tracer for cool molecular 
gas. 

A detailed look at some CO spectra from the 
NGC 6334 FIR II region, presented in Fig. [TO! reveals ve- 
locity components in addition to the main feature dis- 
cussed in Sect. 3.1 and displayed in Fig. [2] Beside this main 
component at -6.5 to -2 km s _1 , there is a weak +6 km s _1 
feature that is only seen in CO and that is present in the 
N, E, S, and center. Apparently, this component is related 
neither to the molecular void nor to the ridge of dense 
molecular gas in the southwestern part of our map. In the 
north, there is a more notable component with a velocity 
of about -lOkrns -1 . It is best seen in the 13 CO(2 — 1) line 
(Fig. IA.lj) . but is also prominent in CO (Fig.[T0]), HCN, 
and HCO + , as well as in [C i] near the northern edge of the 
cross (Fig. [3]). This component marks the southern tip of 
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Fig. 10. CO J=l-0 spectra near the molecular void (for 
the full area, see Fig. IA.lj) showing red- and blueshifted 
emission in addition to the main component. The dis- 
played velocity range is -20kms _1 < Vlsr < +20kms~ 1 . 
The T mb range is -5K < T mb < +60 K. The reference 
position is a 1950 = 17 h 17 m 3O!0, <5i 950 = -35°47'30" 
(a 2000 = 17 h 20 m 51 i ;2, (5 2000 = -35°50'28"). The channel 
spacing is O.llkms -1 . 



a kinematic feature extending southwards from NGC 6334 
FIR I (see Fig. 11 of Kraemer & Jackson [T9"9"9|> . 

There also seems to be some evidence for a molec- 
ular outflow. Redshifted linewings of the 'main' feature 
are apparent near the southwestern edge of the molecu- 
lar void, not only in CO but also in 13 CO, HCN, and 
HCOA Figure fTTI shows the CO map of this redshifted 
emission. The redshifted linewing extends from the stellar 
sources IRS-23 and IRS-24 southwards, along the west- 
ern edge of the molecular cloud ridge (see Fig.[l| that 
shows high (column) densities (Sect. 3.1). The correspond- 
ing blue-shifted linewings are, however, missing (Fig. llip . 
Furthermore, there is no indication of an outflow associ- 
ated with the vibrationally excited H 2 gas observed by 
Straw k Hyland (|1989b|) . 



4.7. ls O/ 17 ratios 

The determination of CNO isotope ratios plays a key role 
in our understanding of stellar nucleosynthesis and chemi- 
cal evolution (e.g., Wilson & Rood ll99~4")) . One of the most 
interesting CNO isotope ratios is that of the two rare oxy- 
gen nuclei, 18 and 17 0. With 18 being mainly synthe- 
sized in massive stars (~8-2OM0) and 17 originating 
predominantly from lower mass stars (e.g., Hoffman et al. 
2001} Stoesz & Herwig |2003p . their abundance ratio may 
become a tracer of material ejected by high-mass stars. 
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17 17 30 25 20 

Fig. 11. Map of the red- and blueshifted CO compo- 
nents. Velocity ranges are -14kms _1 < Vlsr < -8 km s" 1 
for the blue linewings (solid contours) and 2kms~ 1 < 
Vlsr < 12kms _1 for the red linewings (dashed contours). 
Contour levels are 30%-90% of the integrated peak emis- 
sion of 51Kkms -1 (blue) and 41Kkms _1 (red) in steps 
of 10%. The two hexagonal dots mark the stellar sources 
IRS23 and IRS24. 



Observing the J=l-0 lines of C 18 and C 17 has 
the advantage that both lines tend to be optically thin 
and that measured ratios are small (<10) so that re- 
quired observational sensitivities are not extremely dif- 
ferent. Fractionation due to ion exchange reactions (see 
Watson et al. 119761 for the stable carbon isotopes) does 
not play a major role (Langer et al. I1984|) . Hence the inte- 
grated intensity ratio is a good measure of 18 0/ 17 ratios 
(see Wouterloot et al. 2005 for an extensive multi-level 
analysis). While the ratio is small (1.6±0.3) in the Large 
Magellanic Cloud, in the nuclear starbursts of NGC 253 
and NGC 4945 18 0/ 17 ratios of -6.5 have been mea- 
sured (Heikkila, et al. UWEl Harrison et al. \MM Wang et 
al. I2004|) . These results indicate that the ratio does not 
only depend on the relative importance of massive versus 
less massive stars, but that metallicity also plays a role. In 
the Galaxy, the ratio is of order 4 (Table H|), but there is a 
clear discrepancy with the solar system (5.5; e.g., Anders 
& Grevesse[T989) that may have been formed out of mate- 
rial enriched by nearby massive stars (e.g., Lyons & Young 
12005)1 . 

From our spectral data integrated intensities were ob- 
tained by Gaussian fitting or, in cases of poor signal-to- 
noise ratios, by the first moment (area). In the latter 
cases we determined integrated intensities by summing 
line intensities over the velocity interval -13 to +3kms -1 . 
C 18 0/C 17 ratios were calculated for those 152 of the 169 
positions that show C 17 signal-to-noise ratios ^3. 




I(C 0)(K km/s) 



Fig. 12. C 18 0/C 17 integrated line intensity ratios ver- 
sus C 17 J=l-0 integrated intensity. 



In Fig.[TJl the ratios are plotted against the corre- 
sponding integrated intensities of C 17 0. The ratios show 
significant scatter at low integrated intensities due to 
noise in the C 17 spectra. For clouds filling the beam, 
total C 18 and C 17 column densities are proportional 
to v~~ 2 times the integrated line intensity (e.g., Linkc 
et al. 1977). Measured line intensity ratios thus have to 
be corrected: C 18 0/C 17 0(corrected) = (^ci7o/^ci8o) 2 
x C 18 0/C 17 0(observed). The correction factor is 1.047. 
Applying this factor, we obtain abundance ratios of 
4.32±0.13 averaging all plotted data points, 4.13±0.13 
for the positions with 7(C 17 0)>lKkms" 1 , 4.02±0.10 for 
those with 7(C 17 0)> 2Kkms~\ and 4.13±0.12 for those 
with J(C 17 0)>3Kkms _1 , the errors denoting the stan- 
dard deviation of the mean. While the number of values 
obtained toward positions with increasing minimum in- 
tensity is decreasing, the scatter in the individual ratios 
is also decreasing because of higher signal-to-noise ratios. 
Thus standard deviations remain almost the same. 

The derived values are in excellent agreement with 
those of the LTE and LVG (large velocity gradient) model 
results presented by Wouterloot et al. ()2005p from the 
Ophiuchus cloud, with that reported by Bensch et al. 
pUUTj) from a 13 C 18 0/ 13 C 17 intensity ratio toward the 
same cloud complex (see Table 2]), and with those of Ladd 
et al. (2004j) from the similarly nearby Taurus molec- 
ular cloud. Our frequency-corrected values are higher 
than those of the "only" frequency-corrected values of 
Wouterloot et al. (|2005p . but this is not a serious dis- 
agreement. Wouterloot et al. (2005) purposely observed 
positions with a wide range of column densities, including 
some positions where the C 18 opacity was not negligible. 
In these cases, the linear correlation between line strength 
and column density breaks down and abundance ratios 
become larger than frequency-corrected line intensity ra- 
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Table 4. Measured 18 0/ 17 ratios 3 -) 



Source 



Ratio 



Notes 



Reference 



Galactic disk 3.65±0.15 

Nearby dark clouds 3.2±0.2 

Ophiuchus cloud complex (1 position) 4.15±0.52 

Taurus molecular cloud (648 positions) 4.0±0.5 

Ophiuchus cloud complex (20 positions) 3.53±0.11 

Ophiuchus cloud complex (17 positions) 4.07±0.32 

Ophiuchus cloud complex (21 positions) 4.10±0.14 

NGC6334FIRII (169 positions) 4.13±0.13 



frequency corrected value 
frequency corrected value 
frequency corrected value 

from a "simple model" 
frequency corrected value 
assuming LTE 
from an LVG model 
frequency corrected value 



Penzias CEED 
Wilson et al. (fl98T|) 
Bensch et al. ((2001]) 
Ladd flUOU) 
Wouterloot et al. l|2005j> 
Wouterloot et al. ([2005) 
Wouterloot et al. ((2005) 
this paper 



a) Results are mainly derived from C ls O/C 17 J=l-0 integrated intensity ratios. Exceptions are Bensch et al. (|2001f) . who 
measured a 13 C ls O/ 13 C 17 J=l-0 ratio, and Wouterloot et al. PUUS]) . who analyzed the J=l-0, 2-1, and 3-2 lines of C 18 
and C 17 including some lines-of-sight with moderately optically thick C ls O emission. 



tios. This effect may also play a role in the data presented 
by Penzias (fTMT[) and Wilson et al. (fTMTJ) . 

Our data do not hint at optical depth effects. The 
C 18 0/C 17 ratio stays constant up to the positions with 
highest line intensity. 18 O enriched material from mas- 
sive stars and isotope selective photodissociation would 
increase the measured line intensity ratios. High opaci- 
ties in C 18 would reduce them. Only having measured 
the J— 1-0 lines of C 18 and C 17 0, we cannot firmly ex- 
clude that these three effects precisely cancel each other. 
However, in view of our estimate of r( 13 CO J=l-0) ~ 0.3 
in Sect. 4.4 and following the principle of Occam's razor, 
such a scenario is farfetched. Line saturation and isotope 
selective photodissociation can thus be neglected. 

The aim of the isotope ratio measurements was to find 
out, whether 18 enrichment can be seen at the very edge 
of the molecular void. This is neither apparent in Fig. [12] 
nor is there any such tendency when only analyzing those 
positions that are located next to the Hn region. At the 
eastern edge of the void, averaging spectra at positions 
with offsets (0,40), (0,60), (0,80), and (0,100) (see Figs.QI] 
and lA.lj) . we find a frequency-corrected ratio of 3.0±0.2. 
Toward its western edge, at positions (-80,40), (-80,60), 
(-80,80), (-80,100), and (-80,120), we find 2.3±0.1. The 
latter ratio appears to be too low because toward one 
position, C 17 but not C 18 0, is detected. Thus point- 
ing uncertainties might play a role. The ring with radius 
40"-50" and density ra(H 2 ) ~ 10 3 cm~ 3 around the cen- 
ter of the void (representing the inner edge of the cloud, 
see Sects. 3 and 4.4) contains about 5-10 M©, which is 
similar to the mass formerly filling the ionized bubble, if 
we assume the same average density. Depending on stel- 
lar mass and evolutionary stage, the gas at the edge of 
the void might therefore be significantly affected by the 
ejecta (possibly also amounting to several M©) of young- 
massive stars. We note, however, that the 18 0/ 17 ratios 
near the molecular void suggest the presence of lower, not 
higher 18 0/ 17 ratios when comparing them with values 
from interstellar clouds in the solar neighborhood. This 
does not agree with the high values known to characterize 
extragalactic nuclear starbursts. 



4.8. Formation and evolutionary stage of the molecular 
void 

As already indicated, in nuclear starbursts of nearby 
galaxies 18 O is enriched with respect to 17 O (Harrison et 
al. [T999l Wang et al. 12004) . In the Orion OB association 
the youngest sub-associations show higher oxygen abun- 
dances than the older sub-associations (Cunha & Lambert 
1992}. The solar system shows a significantly enhanced 
18 0/ 17 ratio relative to interstellar values, presumably 
because of 18 enrichment just prior to its formation 
(Henkel & Mauersberger 119931) . Why do we not measure 
enhanced 18 0/ 17 ratios near the edge of the NGC 6334- 
FIRII molecular cloud? Following Heger & Langer J2000) , 
Hoffman et al. 1(20010 , and Rauscher et al. (|2002|) . stars 
with zero age main sequence masses up to 20 M© show 
high 18 yields during core collapse. With a highly ob- 
scured 06.5 star (IRS-24; M~40M©) likely being respon- 
sible for the molecular void (Sect. 4.1), stars with ^20 M© 
had, however, not yet sufficient time to reach the super- 
nova stage of evolution. For rotating massive stars, Heger 
& Langer (|2000|) predict enhanced 17 and depleted 18 
surface abundances during the pre-supernova stages, at 
the time of the first dredge up and even before, at core hy- 
drogen exhaustion. Prior to core collapse, a low 18 0/ 17 
ratio in the affected interstellar medium is therefore con- 
sistent with theory. 

The lifetime of an 06.5 star is ~4.5xl0 6 yr (e.g., 
Massey 2005), providing an upper age limit for the molec- 
ular void. There are also indications that the region is not 
extremely young (Sect. 4.1). The extent of the FIR emis- 
sion and the size of the radio continuum source (A16 cm 
diameter: ^lpc; e.g., Brooks & Whiteoak l200ip are large 
when compared with the other main sources of activity in 
NGC 6334. Remarkable is the regular structure of the ra- 
dio continuum shell. Apparently, while being very clumpy 
on small scales (Sect. 4.5), the surrounding parental molec- 
ular cloud must have been quite homogeneous on larger 
scales. There is no hint for the presence of a "champagne 
phase" expected when the ionized shell expands to the 
edge of the cloud. 
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As already indicated in Sect. 1, Hn regions expanding 
into an approximately homogeneous medium should show, 
at later stages of their evolution, a decelerated subsonic 
"D-type" ionization front that is slowly following a shock 
front farther out. This is likely what we are viewing. While 
thermal continuum radiation encompasses an almost cir- 
cular area of diameter 2', even further to the east there is 
the prominent arc of warm presumably shocked H2 gas. 
Since CO is not providing any evidence for outflowing ma- 
terial, H2 may not trace the location of an outflow (as 
suggested by Straw & H viand II 989bp . but instead marks 
the shockfront compressing neutral gas at the periphery 
of the evolved H 11 region surrounding IRS-24. 

In Sects. 4.4 and 4.5 we obtained cloud densities of 
10 3 cm -3 from CO and ^10 5 cm -3 from HCN. For an in- 
termediate value, n(H2) = 10 4 cm -3 , Fig.[4j: provides the 
size of the "initial" Stromgren sphere, 0.35 pc or 0'.7. This 
is smaller than the observed sizes of the molecular void and 
H 11 region, ~1 '.5 and 2' (Figs.Q]and[4|). Since an H 11 region 
can expand up to about five times its Stromgren radius 
(until pressure equilibrium is attained) if the ionizing star 
does not move off the main sequence before (e.g., Yorke 
1986), there is no serious contradiction between models 
and observations. In case the ionized gas just occupies the 
"initial" Stromgren sphere, we obtain an average density 
of n(H 2 ) ~ 2xl0 3 cm" 3 (e.g., Tielens 2005, p. 231). This 
is consistent with our density estimate from CO, which 
represents the bulk of the molecular gas surrounding the 
void. 

In between the outer shell of shock heated vibra- 
tionally excited H2 (Straw & Hyland I1989b|) and the 
ionized region further inside, conditions should be op- 
timal for PDRs. While [C 1] is weak in the molecular 
void that is readily observed in the radio continuum 
through thermal free-free emission, [C 11] is quite promi- 
nent (Kraemer et al. 120000 . Suitable discriminators be- 
tween PDRs and shocks are the absolute and relative in- 
tensities of the [C 11] 158 /zm and [0 1] 63 /im fine structure 
lines. One can find weak [C 11] and strong [0 1] emission 
(7([Oi])//([Cii])>10) in shocks, while the reverse holds 
for PDRs (Genzcl 1991). Observations of far-infrared fine- 
structure lines in NGC 6334, [C 11] 158 /im, and [0 1] 63 /im 
and 145 /im, strongly support the presence of PDRs that 
must be widespread. The most intense [C 11] emission in 
NGC 6334 is associated with the FIR II region (Boreiko & 
BetzHMSl Kraemer et al. [TM51 IMKUl . I([0 l])/7([Cll]) ~ 
3 (Kraemer et al. 1998). 

5. Summary 

Having performed observations of molecular and atomic 
lines in the galactic star-forming region NGC 6334, we find 
the following main results: 

(1) The large-scale distribution of 13 CO J=l-0 emis- 
sion over the entire star forming-region shows strong peaks 
toward FIR I, IV, and V, while FIR II and III are merely 
associated with ridge emission extending southwards from 
FIR I and northeastwards from FIR IV, respectively. 



(2) A detailed study of the eastern part of the most lu- 
minous source of far-infrared emission, NGC 6334 FIR II, 
reveals a molecular void associated with the H11 region 
NGC6334-D. Molecular and [Ci] excitation indicate Tki n 
~40-50K in the surrounding molecular cloud, while to- 
ward the void itself Tki n 200 K. Molecular column den- 
sities and measured obscuration suggest that the exciting 
star IRS-24 (presumably of type 06.5) is located behind 
the bulk of the molecular column. 

(3) The CO spectra show redshifted velocity wings 
south of IRS-24, but the corresponding blueshifted fea- 
ture is missing. 

(4) ls O/ 17 abundance ratios reflecting the 18 0/ 17 
isotope ratio were obtained from measurements of the 
J=l-0 lines of C ls O and C 17 0. The mean ratio is 
4.13±0.13, in good agreement with recent measurements 
toward nearby molecular clouds. At the edge of the molec- 
ular void, however, ratios might be smaller. 

(5) The regular morphology of the ionized gas of the 
H 11 region indicates the presence of a rather homogeneous 
parental molecular cloud on parsec scales. The D-type ion- 
ization front (from observations of the radio continuum) 
and presumably also the shock front farther out in the 
molecular gas (from measurements of vibrationally excited 
H2) can be identified. The age of the Hn region must be 
<5xl0 6 yr- 

(6) The [C 1] distribution follows that of the molecular 
gas and does not show a bright rim at the inner portion of 
the molecular shell. This may be interpreted in terms of a 
high degree of small-scale clumping, with PDRs reaching 
deep into the dense surrounding cloud. Such small scale 
clumping is indicated by the different densities obtained 
from CO and HCN data (n(H 2 )-10 3 and ^10 5 cm" 3 , re- 
spectively) . Observationally, however, the [C 1] lines do not 
provide any direct hint for a close connection with PDRs. 

So far, the velocity field of the shocked H 2 gas has 
not been measured. Additional observations of tracers of 
shocked molecular gas like SiO (Gueth et al. 1998) or 
HNCO (Zinchenko et al. |2"UDU1) would also be desirable. A 
dedicated experiment, investigating the potential anomaly 
in the 18 0/ 17 ratio near the edge of the void to further 
constrain enrichment by ejecta from massive stars would 
be another worthwhile project. 
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Fig. 4. Finder charts of the features on a larger (a, b) and smaller (c) scale toward NGC 6334. The inset in a) (white) 
and b) (black) shows the region sketched in more detail in c). The grey scale outlines the integrated 13 CO J=l— 
intensity (this paper), a) Contours show the distribution of the 1mm continuum emission (Gezari & Blitz 1990). 
Triangles with Roman numbers mark the positions of the 71 /im peaks (Loughran et al. 1986). b) White contours 
show the distribution of shocked H2 gas (Straw & Hvland f~1989bj) . with regions of strongest emission being emphasized 
by white diagonal lines. Prominent Hn regions are marked by A, B, C, D, E, F (see Rodriguez et al. 1982). c) The 
Stromgren sphere of the Hn region D is shown as a black dashed ellipse for n(H2) = 10 4 cm~ 3 (i? s tr oc n(H2)~ 2 / 3 ; see 
Sect. 4.8). The two if-band stars IRS-23 (early B) and IRS-24 (likely 06.5) are given as smaller and larger black 
asterisks, respectively (e.g., Straw et al. 1989c]) . The 1mm peak identified by Gezari & Blitz (|1990[) is plotted as a 
small dot in the southwestern part of the map. The 42/zm FIR peak (Loughran et al. [1.986) is marked as a black square 
with white error bars and the H2O maser identified by Moran & Rodriguez (|1980[) is displayed as a cross, also with 
white error bars, near the northwestern edge of the map. 
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Fig. A.l. Spectra of seven molecular transitions across the NGC6334 FIR II region. Left (from top to bottom): 
12 CO(1-0), 13 CO(2-l), C 17 0, and C 18 0; Right: 13 CO(1-0), HCN(l-O), and HCO+(1-0). The (0,0) position of the 
maps is ai 950 = 17 h 17 m 3O!0 and 5 lg50 = -35°47'30" (a 200 o = 17 h 20 m 51!2 and S 2 aoa = -35°50'28"). The x-axis of 
each individual spectrum denotes the local standard of rest velocity from -20 to +20kms~ 1 . Main beam brightness 
temperature scales are -5 to 60 K (CO), -1 to 25 K ( 13 CO J =1-0), -5 to 45 K ( 13 CO J =2-1), -0.2 to IK (C 17 0), 



